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Abstract. We study temperature and polarization anisotropics of the cosmic microwave 
background (CMB) radiation sourced from primordial cross-bispectra between metric per- 
turbations and vector fields, which are generated from the inflation model where an inflaton 
and a vector field are coupled. In case the vector field survives after the reheating, both the 
primordial scalar and tensor fluctuations can be enhanced by the anisotropic stress composed 
of the vector fields during radiation dominated era. We show that through this enhancement 
the primordial cross-bispectra generate not only CMB bispectra but also CMB power spec- 
tra. In general, we can expect such cross-bispectra produce the non-trivial mode-coupling 
signals between the scalar and tensor fluctuations. However, we explicitly show that such 
mode-coupling signals do not appear in CMB power spectra. Through the numerical anal- 
ysis of the CMB scalar-mode power spectra, we find that although signals from these cross- 
bispectra are smaller than primary non-electromagnetic ones, these have some characteristic 
features such as negative auto-correlations of the temperature and polarization modes, re- 
spectively. On the other hand, signals from tensor modes are almost comparable to primary 
non-electromagnetic ones and hence the shape of observed S-mode spectrum may deviate 
from the prediction in the non-electromagnetic case. The above imprints may help us to 
judge the existence of the coupling between the scalar and vector fields in the early Universe. 
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1 Introduction 

Owing to accurate observations and data treatments, it has turned out that galaxies and 
clusters of galaxies hold micro-gauss magnetic fields (e.g., [1-4]). Furthermore, studies on 
some astrophysical processes suggest that the strength of magnetic field in the inter-galactic 
medium is larger than O(10~^^ — 10"^'') gauss [5-8]. What is the origin of these large-scale 
magnetic fields? It may be generated via an interaction between infiatons and some sort 
of vector fields, which breaks the conformal invariance, in the inflationary Universe (e.g., 
[9-13]) ^ and such kind of model can generate nano-Gauss magnetic fields at the present 
Universe [12], which is comparable to upper limits from the power spectra [18-29], bispectra 
[30-37] and trispectra [38] of CMB anisotropies. However, the models with such coupling 
might suffer from some problems. One is that the energy density of the vector field spoils 
the inflationary background dynamics, so-called, backreaction problem, and another is that 
due to the time-dependence of the coupling the Universe enters strong coupling regime at the 
early stage of the inflation, so-called strong coupling problem (e.g., see ref. [13]). Hence, it 
has been commonly understood that constructing a realistic model of generating primordial 
magnetic fields during infiation is quite difficult. There are also several works about universal 

^ Other than the inflationary origins, diverse generation and amphfication mechanisms of the magnetic field 
have been proposed (e.g., [14-17]). 
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bounds for the inflationary model from the current observations of the large scale magnetic 



Recently, such vector field and coupling have been also widely discussed beyond the 
context of magnetogenesis (e.g., [42-45]). Because, the interaction between the inflaton and 
the vector field can induce the cosmological three-point cross-correlation (cross-bispectrum) 
between the fluctuations of the scalar field, metric perturbations and electromagnetic fields 
and there must be some cosmological imprints of it. 

In ref. [46, 47], although in the context of the generation of the magnetic field, the au- 
thors computed such a cross-bispectrum composed of one scalar field and two magnetic fields, 
and analyzed its shape dependence in the similar manner as the cases of non-Gaussianities 
without magnetic fields [48-50] . They suggested that this cross-bispectrum may be accessible 
through the cosmic microwave background (CMB) observations and a combined survey of 
large scale structure and Faraday rotation [51]. Therefore closely estimating CMB signals gen- 
erated from the cross-bispectrum will be an interesting and important work for studying the 
early Universe. Furthermore, although refs. [46, 47] have not discussed, the cross-bispectra 
involving tensor perturbation or the electric field can also be generated and hence may cre- 
ate characteristic CMB fluctuations. In refs. [45, 52], the authors investigate the primordial 
fluctuations sourced from the vector field through the coupling during infiation and estimate 
the non-Gaussianity of the primordial curvature fluctuations. 

This paper examines impacts of such electromagnetic-scalar and electromagnetic-tensor 
bispectra on CMB anisotropics. Electromagnetic fields, which are created during infiation, 
induce the anisotropic stress fiuctuation (depending quadratically on electromagnetic fields) 
and affect the subsequent evolution of the CMB fluctuation of scalar, vector and tensor modes 
even during the radiation dominated era [19, 20, 24]. This implies that the electromagnetic- 
scalar and electromagnetic-tensor bispectra become sources of not only CMB bispectra but 
also CMB power spectra. Hence, flrst, we construct complete formulae for these cross- 
bispectra and the power spectra of temperature and polarization anisotropics. In such case, 
since the scalar and tensor modes sourced from such anisotropic stress fluctuations have 
same origins, it could be readily imagined that the scalar-tensor mode coupling appears. 
However, we explicitly show that the CMB power spectra do not have such mode-coupling 
components. Then, we find that CMB power spectra from the electromagnetic-scalar and 
electromagnetic-tensor bispectra have characteristic features. In the power spectrum of the 
scalar modes, although these are subdominant contributions compared with the primary non- 
electromagnetic signals [50] , the auto-correlations of the temperature (intensity) and E-mode 
polarization have negative values, respectively. CMB power spectra induced from the tensor 
modes have almost the same amplitude of primary ones. In this sense, the effects on the 
-B-mode spectrum might be quite interesting. 

This paper is organized as follows. In the next section, we summarize some behaviors 
of the vector field and metric perturbations in the context of slow-roll infiation. In section 3, 
we present the formulae for the electromagnetic-scalar and electromagnetic-tensor bispectra. 
In section 4, we compute and analyze the CMB signals generated from these cross-bispectra. 
The final section is devoted to the summary and discussion of this paper. 

Throughout this paper, we obey the definition of the Fourier transformation as 



and the rule that the subscripts and superscripts of the Greek characters and alphabets run 



fields [39-41]. 
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from to 3 and from 1 to 3, respectively. 



2 Nature of inflation coupled with the vector field 

In this section, let us discuss some background and perturbative quantities in a simple in- 
flationary model with the coupling between a scalar field (inflaton) and a vector field. The 
action is given by 



S = (Tx 



Ml, 1 



+ Sv , (2.1) 



Sv = - j d^x^^g''^g''''W{ip)F^,Fx. , (2.2) 



where (f is the inflaton, F^jy = d^A^ — dyA^ with A^ being the vector field, and Mpi = l/yG 
is the Planck mass. For the dynamics of the vector field, only the time-dependence of the 
coupling function is important. Here, we simply impose the power-law type as 

W^(r/) = W^/(^)'" , (2.3) 

where a{rj) is the scale factor with r/ being the conformal time, and the subscript / means 
the value at the end of inflation, rj = r]j. This power-law shape is often realized by an 
assumption of a dilaton-like coupling [47] . Due to the time dependence of this coupling, the 
conformal invariance of the action (2.2) is broken and hence the vector field can be survive 
even in the inflationary expansion. In addition, let us consider the case that after the end 
of inflation, W = Wj is kept and the conformal invariance is restored. This seems to be 
natural assumption, since the coupling function W is, here, dependent on the inflaton and in 
general the inflaton falls into a stable state after the end of inflation. Owing to this condition, 
we do not need to consider the additional growth of the vector field at late stages, which 
is disfavored by observations. In the literature, to identify such the vector field with the 
observed magnetic field, the coupling function at the end of inflation Wj should be taken to 
be unity [9-13]. 

The evolution of the vector field strongly depends on the power-law index of the coupling 
function, n. In order not to spoil inflation due to the backreaction of the energy density of 
the vector field, we have a bound on the spectral tilt of the running coupling: |?7,| < 2.1 (e.g., 
see [12, 13, 47]). 

This type of magnetogenesis model has another problem. In case where n > and Wj 
is taken to be unity, the running coupling constant, W{r]), is so small at the beginning of 
inflation and it means that the gauge coupling at the beginning of inflation is quite large, 
that is, inflation starts at the strong coupling regime. In such regime, we can not treat the 
perturbation of the vector field as a free field any more. However, refs. [46, 47] pointed out 
the possibility to avoid such a strong coupling problem by the action arising from the UV 
completion and including the violation of the gauge invariance as described in theories with 
extra dimensions. On the other hand, in ref. [45], the authors mentioned another possibility 
that the vector field localizes in a hidden sector and generates real electromagnetic field 
through any hidden coupling. In the following discussions, although the solution of the strong 
coupling remains unclear, we shall admit the above treatments and continue to discuss about 
cosmological imprints of the vector field. 
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In this paper, we assume the single field slow-roll inflation model where the slow-roll 
parameters, e = —H/H^ and 5 = (p/{Htp), are almost constant and much less than unity. 
Here H is the Hubble parameter and ' = d/dt denotes the derivative with respect to the 
physical time. Then, we have the relations between background quantities as 

aH 1 — e aj \rji J 
2.1 Evolutions of the fluctuations 

Let us evaluate the evolution of curvature and tensor perturbations, and the fluctuation of the 
vector field, following ref. [47]. Using the Arnowitt-Deser-Misner formalism, each component 
of the scalar- and tensor-mode metric up to the first order is given by ^ 

goo = -N^ + QijN'W , go^ = QijN^ , g^J = a^[il + 2n)5ij + h,,] , (2.5) 

with 

Ar=i + ^, N^ = -^+ed,V-^n. (2.6) 
H aH 

Here, we obey the gauge, 6ip = 0, and the transverse-traceless condition, d^hij = = 0. 
These metric perturbations and the vector field are quantized as 



7^(x,r/) = / ^ [7^fc(r?)e^l^-/3(k,0) +h.c." 
M^,V)= /tC^ E [^fc(^)e^''"0f^(k)/3(k,A) + /i.c.] , ^2.7) 

' A=±l 

where /?(k. A) and /3^(k, A) are respectively the annihilation and creation operators of 7^ (A = 
0), Ai (A = ±1) and hij (A = ±2), and satisfy /3(k,A)|0) = and [/3(k, A), /?t(k'. A')] = 
{27r)^6x^x'd{k — 'k'). In terms of the vector field, we have adopted the Coulomb gauge d'''Ai = 

and Ao = 0. The explicit forms of the divergenceless vector, O^-^^^ and transverse-traceless 

(±2) 

tensor, O-^ , are shown in appendix A. Then, as the mode function of each perturbation, 
we have [47, 53] 

TZkiv) = 'Jlk*Uu{-kr]) , 
Vk{l]) = Vk*Ua{-kl]) , (2.8) 

where we have defined the time-independent and time-dependent functions, respectively, as 
TZk. = ie-'^/2+-/4(l - e)2^m-^^tll!l^^ , (2.9) 



For decaying nature, we neglect the vector-mode metric perturbation. 
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= _2ie-'^/2+-/4(i _ ,)2Mr(/.)^^-^^^^ , (2.11) 

and 

^-(^) = ^^^^^'^(^)' (2-12) 

with T(m) and Hm\x) being the Gamma function and Hankel function of the first kind, 
respectively. Here, the subscript * denotes the value at a given time after the scale of interest 
exits the horizon, = With the asymptotic form of Um{x) at —kr]^ 0, we find that on 
superhorizon scales TZk{r]) ~ TZk* and /ia:(??) — ^fc*- Properties of this function are described 
in appendix B. Subscripts of this function are determined by the slow-roll parameters and 
the spectral index of running coupling as 

3 e In e + (5 

^=2 + 137' "=2 + T^' ^ = ^^+137- (2-13) 

2.2 Power spectrum of each mode 

By use of the above equations, the power spectra of these perturbative quantities are imme- 
diately calculated: 

(7^(k,7?)7^(k^7?')> = (27^)37^,(,?)7^*(r?')5(k + k') , 
{A{k,r^)A,{k',r^')) = i2nfvk{r])vUr]')Pij{mk + k') , (2.14) 
(hijik,7])hMik',7j')) = i27rfhki7])hUr^')Uij,M6ik + k') , 



with 



P,,(k)^ e]'\k)ey'\k) = 6ij - kk 

A=±l 

n,,,«(k)= Ye^{k)ei-'\k). 



(2.15) 



A=±2 



One can find conventions of the polarization vector and tensor, e^^^ and e[j^ in appendix A. 
The vector field can be decomposed into the electric and magnetic components as 

Ei(k,77) = -^Kk,r?) , ^^^^^ 
Bi{k,r]) = ir]ijkkjAkik,r]) , 

where ' = d/drj denotes the derivative with respect to the conformal time and rjabc is the 
3D antisymmetric tensor normalized as r/123 = 1. Hence, the electric and magnetic power 
spectra outside the horizon at the end of inflation, namely, —kr]j <C 1, are respectively given 

by 

{Ei{k,vi)E^{k',vi)) = (2vr)3:^P,,(k)5(k + k') , 
{B,ik,r^j)B^ik',vj)) = (2vr)3^^P,,.(k)5(k + k') , 
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The magnetic component outside the horizon at the end of inflation re-enters the horizon 
at late stages and may become observed magnetic fields in case the vector fields can be 
considered as the direct seeds of the present magnetic fields. On the other hand, in such 
case the electric component quickly decays due to the large conductivity (see, e.g., [12]). 
According to ref. [47], with the WMAP results [50] and the assumption that e = 0.01 and 
the instantaneous reheating, the strength of the magnetic field at the present time can be 
calculated as 



^ , 2v ^ k^PBik)aj 
dlnk^ ^' 27r2 

^ ^0-22.8-22.5nB ^ I 2 ) ( ^ \ QaUSS^, (2.18) 

WiT%^ VMpc-V 



where ub = 4 — 2n/(l — e). This equation expresses that the magnetic field is drastically 
suppressed as ns increases. Hence, we can see that for n = 2.1, i.e., ns = —0.2, and Wj = 1, 
the strongest magnetic field is produced as -BiMpc — 10~^Gauss. This value is comparable 
to upper bounds obtained from current CMB observations [23, 26, 27, 29, 34-38] and hence 
we can expect that the electromagnetic-scalar and electromagnetic-tensor correlations also 
have detectable impacts on CMB fiuctuations. Thus, in section 4, we focus on the analysis 
of CMB signals for ns = —0.2, i.e., n = 2.1 and in order to treat more generic case we also 
study the effect of the electric component without considering the dilution of it. 



3 Primordial cross-bispectra 



In this section, we calculate the primordial cross-bispectra between two electromagnetic fields 
and one curvature perturbation {{TZEE) and (TZBB)), and those between two electromag- 
netic fields and one gravitational wave {{hEE) and {hBB)) ^. 

As a powerful tool to compute the higher-order correlation in the inflationary Universe, 
we utilize the in-in formalism [54, 55]. In this formalism, the expectation value of the time- 
dependent operator in the interaction picture, 0{t), is given by 



m)) = (o 



(3.1) 



where T and i/int(t) are the time-ordering operator and interaction Hamiltonian, respectively. 
To compute the tree-level bispectra of primordial fluctuations, it is only necessary to use the 
first-order expression as 



(0(7?)) = -2Im r dr]' (0 | : Fint(?/) : 0(r?) | 0> , (3.2) 

J — oo 

where : : denotes the normal product. In the following discussion, we equate {TZEE), 
{TZBB), (hEE) and (hBB) to {0(1])). 



^We do not deal with (TZEB) and (hEB) because these correlations do not generate the CMB signals as 
seen in section 4. 
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3.1 Electromagnetic-scalar bispectra; {IZEE) and {IZBB) 



Expanding the action (2.2) to linear order in TZ leads to 



S- 



dr]d^xW{if) 



A'/ 



^4 - -^A'jPij + e5i(V"^7e')4F, 



.(3.3) 



Hence, the interaction Hamiltonian of the scalar part is derived as [47] 



n 

n=l 



(2vr)3j^k^ 



\n=l 



7^'(^,r/l) 



1, k',i{6^^'^' - 5^^5^^)A^ik2',m)M^3',m) + (2 o 3)1 1 



a{rji)H{rji 



Here, we have neglected the terms of order e. Substituting equation (3.4) into equation (3.2) 
and overcoming somewhat complicated calculation, one obtain the explicit form of the 
electromagnetic-scalar bispectrum {Y = E, B) as 



(7^(kl,^?)r(k2,r?)y,■(k3,?7)> = K^k^ 



\n=l 



,n=l 



(3.5) 



where 



K3*l 



Wi{-kir]i 



2n 
1-e 



(3.6) 
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A'^^^'s and A/^^'s involve the time integrals and angular dependence arising from contrac- 
tions, which are respectively given by 

Kf^^ikiMM = -^h Im ul{-hr])u':{-k2r])u':{-ksv) 

r-CO _2r^ / d \ d d 

X j dTiT^ ( 1 + (1 - e)ri— J u,,(ri)— nQ(j;2n)^itQ(a;3Ti) , 

Kf^^{k^MM = 2^ Im n:(-A:ir/)<*(-A;2??X(-fe3^) 

f°° / d \ 

X J dTiT^ {l + {I - e)Tl — j UuiTi)Ua{x2Ti)Uaix3Ti) , 

Kf^'^ik^MM) = 2^ Im ul{-k,r^)u'*{-k2il)u':{-k3r^) 

X / drir^ ' ' (1 - e)u;.(ri)ua(x2ri)— nQ,(x3Ti) , 

K^''''{ki,k2,k3) = K^^^{ki,ks,k2) , 

Kf-^'^'iki, k2, ks) = -Akik2ks Im n:(-A:i7?X(-A:2r?)n* (-fcsr?) ^^'^^ 

driT^ ' M 1 + (1 - e)Ti— Mz.(ri)— tia(2;2Ti)— ^^(xsri) , 
kiT] V ""^1/ "''"1 "■'"1 

K^''''{ki,k2,ks) = 2MM Im <(-A:ir?X(-A:2?7)<(-A:3r/) 
ki 

f°° ( d \ 

X y driT^ ' ' [^1 + (1 - e)Ti— J M,.(ri)iiQ,(x2ri)tiQ,(x3Ti) , 

ii'3^''''(A:i,A;2,A:3) = 2^^ Im <(-A;ir?)<(-A:2?7)<(-A:3r?) 

X / driT^ ' '(1 - e)'Ui/(n)ifa(a:;2Ti)— tiaCxsTi) , 

A:f^^(A;i,A;2,fc3) = i^3'^''''(fci,fe3,A;2) , 
and 

A/;^(2) ^ . k3P.fc(k2)Pjfc(k3) + fc2/4fc^'^fc(k2)P,7(k3) , 

A/;f' ^ = - + k2 • ks) Pa.(k2)P,fc(k3) + 4^fc3fcP^fc(k2)P,7(k3) , 
A/;f' ^ = - + k2 • ks) Pa.(k2)P,fe(k3) + 4^fc3fcP^fc(k2)P,7(k3) , 



A/;f = \ (k2 • k35,, - k2,ku ) , 



(3.8) 



7V^|(^) = (k2 • k3)^(5ij - k2 • k3fc2j4i + VimnhmhnVjklhkksi , 

7V^|(^) = ^k2 • k3 + (k2 ■ ksSij - k2jh?j + mmnkzmk2n'njklhkk2l , 

A/^^(^) = ^k2 • k3 + (^k2 • k35ij - k2jh?j + riimnk2mk3n11jklk2kk3l ■ 
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Here, ti = —kir]i,X2 = k2/ki and X3 = k^/ki. Note that unlike ref. [47], in equation (3.5), we 
do not perform the contraction between electromagnetic fields. It is because off-diagonal com- 
ponents of the electromagnetic-scalar bispectrum create significant signals on CMB anisotropics. 

To see behaviors of these cross-bispectra at the end of inflation on superhorizon scales 
{—kr]i <^ 1), let us focus on two specific cases: n = ±2, i.e., a = 5/2 and —3/2. Then, 
explicit forms of K^^^^s are given in appendix C. Taking into account the r//-dependence of 
K^^^'s and ^(^), we can find a fact: 



{nEE)^_, oc i-krji)-^ , {nBB)^_, « (-kvi)-' 



(3.9) 



This implies that the electric (magnetic) contribution dominates over the cross-bispectrum 
for n = —2 (n = 2). Such inverted behavior between (JZEE) and {TZBB) under the sign 
reversal of n may be observed for any n. 



3.2 Electromagnetic-tensor bispectra; {hEE) and {hBB) 

In the same manner as the scalar case, expanding the action (2.2) to linear order in hij as 



S- 



(h) 



I 



dr]d X 



W=A[A'. - h'^^'idiAkdjAi + dkAAAj) + 2h'^'''diAkdiA 



kl; 



(3.10) 



and transforming into the Fourier components, we gain the tensor-part interaction Hamilto- 
nian: 



:{/i^^(k'i,r?i)A:(k2',r?i)4(k;„77i) 



(3.11) 



+ 



Substituting this Hamiltonian into equation (3.2), we have 



/i,,(ki,77)y*^(k2,r?)ll(k3,r?)) = {2t:)H [Yl^^-j A^\k^M,h:r]) 



where 



ijkl 



.n=l 



(3.12) 



2*1 I'^fea*! 



a4(r?) 



Wi{-kir]i)^-^ 



(3.13) 
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In the tensor case, the time integrals and angular-dependent parts are summarized as 
Kl'^^ = 2h Im u;{-kiri)u':{-k2v)u':i-hri) 

X / dTi ' " U^{Ti)-—Ua{x2Tl)^Ua{x3Ti) , 
J-kif) CiTl OTi 



and 



'■CO 2n 

-kiri 



At^^ = 2A:ifc2A:3 Im u;{-kiri)uU-k2vXi-hv) 

X / dTi ' ' U^{Ti)—Ua[X2Tl)——Ua{X'iTi) , 

/OO 2n 
dri U^,{Ti)Ua{x2Ti)Ua{X'iTi) , 

•^3" = -n„„..j(ki)ft„(k2)fl„(k3) , 

— ^2mfe3nnmn,ij(kl)Pfcr(k2)-Pzr(k3) — ^2 ' ^S^mn,iji}^l)Pkm{^2)Pln{}^S 
(ki) \ PkT{^2)Pln{^z)k2nikzr + Pkn{^2)PlT{^^)k'inik2 

X [VksgVltrhmhn + msmVltnhqhr ^ksq'Hltn 



(3.14) 



(3.15) 



In the same manner as the scalar case, we analyze the cases for n = ±2. From the 
combination of analytic expressions of K^^^^s described in appendix C and A^'^\ we find 
the scaling relations as 



{hEE)^_^ oc (A:r?,)-4 , {hBB),^_^ oc (km) 



4 /.oos n ^-2 (3-16) 



This dependence is identical to the scalar counterpart; hence for positive (negative) n, mag- 
netic (electric) part dominates over the cross-bispectrum. 



4 CMB power spectra 

In this section, we discuss impacts of the electromagnetic-scalar (3.5) and electromagnetic- 
tensor (3.12) bispectra on CMB fluctuations. Then, notice that we do not deal with CMB 
bispectra but CMB power spectra because the CMB fluctuation arise quadratically from 
electromagnetic fields. 
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4.1 Formulation 

CMB anisotropies are quantified by one intensity {X = X) and two linear polarization B) 
fields for the scalar [Z = S), vector (V) and tensor (T) modes ^. These are expanded by the 
spherical harmonics as 

Then, each coefficient is expressed as [56, 57] 

„2 



(4.2) 



where A denotes the helicity of each perturbation: A = (Z = 5*), ±1 iV) and ±2 (T), 
and X discriminates the parity of each field: x = ^ {X =X,E) and 1 {B)^ respectively. ^'^^^ 
expresses the initial perturbation as ^'-"^ = TZ and = = ^O^J'^^htj, and 7^^^ is the 
transfer function derived from the line-of-sight integral. 

As seen in equation (4.2), CMB anisotropies depend strongly on the magnitude of the 
initial fluctuations. Electromagnetic parts of the vector fleld also affect CMB anisotropies 
via the primordial perturbations as follows and those can generate both the scalar and the 
tensor modes. If electromagnetic fields exist in the radiation-dominated era, their anisotropic 
stresses act as sources of scalar and tensor metric perturbations. Due to this, on superhorizon 
scales logarithmically-growing metric perturbations arise prior to neutrino decoupling. After 
this, however, neutrino anisotropic stresses emerge and compensate for electromagnetic ones; 
therefore the enhancement of metric perturbations ceases. Consequently, we have [24, 37] 



AO) 



(k)«i?,in ^o(f (k)n^.,(k) , 

ei±')(k) « 6R,ln iolf )(k)n^.,(k) , 



(4.3) 



where rj^, is the conformal time of neutrino decoupling, ~ 0.6 is the ratio of the en- 
ergy density between photons and all relativistic particles, and a subscript A denotes the 
quantity originated from electric and magnetic parts of the vector field. IlAij means the 
time-independent energy momentum tensor of the residual electromagnetic field after the 
inflation end as 



= E Y^{^',r,j)Y,{\.-).',r,j) , (4.4) 

with and being the photon energy density and its present value, respectively. Note 

—1/2 

that due to oc Wj , H^^. is independent of Wj. This induces the absence of Wj in 



Here, we neglect the circular polarization. 
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The metric perturbations outside the horizon behave as initial conditions of the CMB 
anisotropics of the scalar and tensor modes. Strictly speaking, electromagnetic fields also 
urge the modification of the transfer function. However, this change is negligible at large 
scales [24] and hence we use the transfer functions without depending on electromagnetic 
fields [57-59] in our numerical calculation. 

On the other hand, electromagnetic fields also generate the CMB fluctuation of the 
vector mode because the vector-mode anisotropic stress equates to the Lorentz force and this 
supports the growth of the vorticity at recombination. Thus, the electromagnetic vector mode 
produces characteristic transfer function [19, 20, 35]. As the initial vector-mode perturbation, 
we adopt the fluctuation of the electromagnetic anisotropic stress as [37, 57, 60] 



^off )(k)n^.,(k) . 



(4.5) 



Equations (4.3), (4.4) and (4.5) imply the quadratic dependence of the initial perturbations 
on electromagnetic fields. This means that {IZEE) , (TZBB) , {hEE) and {hBB), which are 
given by equations (3.5) and (3.12), equate to the power spectra of the initial perturbations 
as ^'C^'^^^C^^''^ aiid therefore become sources of CMB power spectra. 

From here, let us focus on a formulation of the CMB power spectra, which is expressed 

as 



\n=l 



(Zu) 



n=l 



(2^)2 



(A2) 



(4.6) 



At first, we should reduce the initial angular power spectra obtained from equations (3.5), 
(3.12), (4.3), (4.4) and (4.5) as 

Wiaj 



c" 



A2 



A-Kp. 



'7,0 



{-If- / d^ki„,.,y,*„jki);,,y,;„jki) 



A^2 d^2, 



X5(ki+k'2 + k'3).^^^^ 



-k2) 



dV3FA,A2(kl,k'2,k'3) 



kl 



(4.7) 



with 

Fo,A,(kl,k'2,k'3 



A^'^\k,,k'„k'„r,j) 

X E i:^n^''''(^i'^2,4)oJ^^^(ki)A/;^'^"^(k^,k;,) 



Y=E,Bn=l 

F±2,A2 (kl, k'2, k'3) ^ A^'^Hkuk'^, k'^,7]i) 

2 



C 



A2 



Y=E,Bn=l 

0) 

(A2 = ±l) . 

[m.,\n{i^J'ni) (A2 = ±2) 



X E E^n''''(^i'^2,fc^)^eS7^^^(ki)Or^(ki)A/;;r(ki,k',,k^^^ 



l(A2). 



.(4.8) 



^ \R^\n{'qy/r]i) (A2 

1 
2 
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Here, we have decomposed the delta function mto 



5(ki+k2) = ^^^i^5(ki + k2). 



(4.9) 



For performing the angular integrals in the Fourier space, it is convenient to expand all 
angular dependent parts by the spin spherical harmonics. Then, the delta function is given 

by 



5(ki 



+ k'2+4) = 8/ y^dy V (-1) ' i ' 3LAkiy)jLMy)jL,{k'^y) 



LIL2L3 
MIM2M3 



■■yi.MA^i)Yl2MAK)Yl,M.M^^^^^ > (4-10) 



where jiix) is the spherical Bessel function and the / symbol is defined by 



JS1S2S3 



(2/i + l)(2/2 + l)(2/3 + l) / h h h 
47r V si S2 S3 



(4.11) 



F\^\^^s also involve the angular dependence. Taking the contractions by use of the conven- 
tions shown in appendix A and refs. [57, 61], we have: 



L'L" 



'=E,Bn=l 

MM'M" 

x(-l) 
where s = A2 — Ai and 

= { 



2 L' L" 



M M' M 



I T\/rii ' 



LL'L 
s+L+L'+L" 



Y=E,Bn=l 



(4.12) 



^2 2 L 



y^ -sYlm (Mi)Yl'M' 0^'2)Yl"M" (k's 



MM'M" 



L L' L" 



^'\M M' M 



I l\,Tll ' 



-2 (A2 = 0) 
2^/3A2 (A2 = ±1) 
^2^3 (A2 = ±2) 



(4.13) 



New functions for Ai = and ±2 are respectively given by 



1 



^UYY _ __ 

^ 2 ^ 



K 



k. 



■K. 



k. 



ko 



K 



Kfyy = -K^yy _ x^yy _ xf^y 
jcf^y = -jcf^y , 



(4.14) 

(4.15) 
(4.16) 



-4/2V(f) ^^^n^/^^aM^^^o + ^L^2)(5L^o + ^L^2){Jf , (4.i7) 
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_ 4 7-01-1 /^^^V .01-1 rOl-1 rOOO rOOO j 2 P2 P3\ j 2 L' L"\ 

^L'L" - "^hll I Y I 2^ -'pall -'psll hp2L'hp3L" 1 | 1 1 Pa P2 ) 



2 L' L" 



^^(3) _ .^01-1 /^47r\ ^-^ ^-01-1 rOl-1 rOOO rOOO J 1 1 I ^4 10^ 

^ ^ P2,P3=0,2 I 1 1 P3 



,;i?(i) rOi-1 /^47ry 



(4.20) 



yB{2) _ 01-1 ( .000 .000 f 2 L" 1 

•^L'L" - "-'211 I Y/ ll-f- 11^" 1 1 111' l^-^^i 



3 



2 L' L" 



and 



•^LL'L" 



V-B{2) 
^LL'L" 



x(5L',o + '5L',2)(5L",o + 5L",2)/°}n ^L''nM 2 1 1 J> , (4.23) 



2 1 1 
2 1 1 



/ Til 



S,/^ /'rOl-l'l^ ^ rOl-1 rOl-1 / 2 P2 P3 \ rOOO rOOO J 9„ „ 

»V3 (^211 ) ( Y / ^ P2ll ^P3ll 1 1 1 1 f hp2L'hp3L" \ 2 P2 P3 

^ ^ P2,P3=0,2 ^ [ill 

+8V3(l--)^(f ) 5: «r/C-/ 2 11 

V / p2,p3=0,2 [ 2 1 1 J 



„ 7-000 7-000 / i\L'+L" 
X^1P,L'^1P3L"(-1) \ 1 p3 p2 



L L' L" 



3 / 

P2,P3 



X [1 + (_i)P.+P3] I 2 P2 } I 2 p2 P3 I , (4.24) 
Vfl^i,, = -2^3 {Ill^'f (^^y 6 I 2 1 1 I 5^, i<^^,, 1 , (4.25) 



-2^^(/^i\"^)'(y) /iTAT'I 2 1 1 



+2v/3 (/°Jr^)^ (^) ' 6(-l)^/??o,/??°„ I 2 I j I { t ^' Y } 
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+2V3 (/°Jr^)^ (f ) ' 4 (^,, + AV^/iT' 



L Lp 1 



Lr 



L Lp 1 



Lp=l,2 



(4.26) 



For calculating the tensor mode, we have used equation (A. 5). By the Wigner symbols, the 
integrals of these spin spherical harmonics reduce to 



■J T III !<, rill 



r-A2A20 
£2 2 V" 



V»M" 



X(-1) 



M" 



I 



h Li L'" \(h2 V" 
mi Ml M'" J \m2 M -M" 



Ai £imi A2 £2m2 -LM LiMi 



ErAlO-Al r 



-A2SA1 
'2 L L'" 



L"'M"' 



X(-1) 



(4.27) 



A/" 



£1 Li L'" 
mi Ml M'" 



£2 L V" 

7712 M -M'" ' ' 



(f Yl^ ¥[, j^j, 



-1)*^ 6l2,l'Sm2 



Furthermore, the summations of Wigner symbols over azimuthal quantum numbers result in 

\M'+M" 



-M' 



E E 

M'M" M2M3 



E (-1) 

M"'MM, 



L V L" \ ( Li L2 L3 
M M' M" ) I Ml M2 M3 



'JL2,L'^M2-M'^L3,L"^M3-M" 



1 



2L + 1 



- 0l2,L'0L3,L" 



(-1)*'^(^Li,l'5mi -m(-I)'^ 



A/+A/" 



ii Li L" 



mi Ml M'" J \m2 M -M 



L L" 



(4.28) 



Sli,L^Mi-M 



(-1) 



-<5Li,L(-l)™'^%/2'^i 



mi ,— m2 



2£i + l 

Here, the scalar- mode counterpart corresponds to the case where L = 2. By the selection 
rules of Wigner symbols [57, 61], we can simply the summations over L'" for Ai = and ±2, 
respectively, as 

5(2£i + 1) 



ErO T- 
^£l2L"'^£l 



A2A2O 
2L'" 



47r 



-'^A2 



T-A1A2-S \^(_^\L"' T-AlO-Al r-A2sAl _ e r-AlAlO/i N^i+Li ST- f T 

-'2 2 Li Z^l ^liLiU'i^hL^L'" — "Ai,A2-'2 2 Lj 

L'" L'" 



AiO-Ai 

eiLiL'" 



(4.29) 
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These results ensure that the couphngs between different circular modes, such as ^^^''''^^"^'^^^^ 
vanish. Due to this fact, the CMB power spectrum of the vector mode does not arise from 
the electromagnetic-scalar and electromagnetic-tensor bispectra. From the above treatments, 
the initial angular power spectra for Ai = and ±2 (4.7) can respectively reduce to 



/2'ni2 

{k2) 



(A2) 



,i2in2 

ik2) 



(-l)-M,„,,5^,_„,5A„o(-3)i?,ln(^^) ("4^) 



X / y^dyj2{kiy) 



vr kl 



2L2L3 



L2L3 



n=2 



n / k'^dk'jLAKy) 



WiA^'^Hki,k'^,k'„in) 



Y=E,Bn=l 



(4.30) 



X E ^2^t 
Li=0,2,4 

00 



y^iiLiL'" 

L»> 

5{ki - k2) 



(2Li + l)(2^i + l) 



xSy^ y^dyiL^{kiy) 



L2L3 



'^i+'^i+Ls 



,n=2 



n / k'^dk'jLAKy) 



WiA^''\ki,k'^,k',,ijj) 



J2 j2K:iyyik„k',,k',)vZl 



Y=E,Bn=l 



(4.31) 



Note that {-l)'^^5t^, 

,£2^mi,—m2 enforces the rotational invariance of the initial and CMB 

power spectrum. 

In order to convert into CMB power spectra, we substitute these equations into equa- 
tion (4.6). As a result, we obtain final formulae for CMB power spectra of the scalar and 
tensor modes, respectively, as 



n (-^n) \ _ ^(^1-22) / -l\mix X 



~f{ZiZ2) 



\n=l 



c 



(SZa) 

XIX2/1 



ry'dy^ I kjdkij2{k,y)Tj,%{k^)Tj,%^iki) 

T T O 



(4.32) 
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^xfxJe-, = 5t,Za ('5xi,o'5x2,o + ^xi,i42,i)6V3i?^ In 



:WiA^^\ki,k'^,k'^,Vi) 



5: E^"^^(^^'^-fc^)vrs 



Y=E,Bn=l 



, (4.33) 



Li=0,2,4 



LL" 



r]ij\ 47r/)^,o 
1 

(2Li + l)(2^i + l) 



x8 / y'dy I kidkijLAkiy)T^x^S^i)T'xV/Ski) 



-(T) 



2^ l-J^j 2 ^LlL2i3 



L2LS 



n 

.71 = 2 



Y^Kfyik,,k',,k',)vlS.„^ 



Y=E,Bn=l 



(4.34) 



In the derivation of equation (4.34), we have performed the summation over Ai = ±2 and A2 
as 



n 



c T-AiAiO / rAiO-Ai 
OAi,A2-'2 2 Li l-'^iLiL'" 



, n=l A„ 



97—220 ( t2 -2 
^-'22Li l^-'^iLiL'" 





(xi + X2 + = even) 
(xi + X2 + = odd) 



(4.35) 



From these expressions, we can see that due to ^(^''^) oc Wf^, CMB power spectra are 
independent of Wj. Therefore, the behaviors and amphtudes of CMB power spectra are 
determined by only the spectral tilt of the running coupling, n, except some inflationary 
parameters. 



4.2 Numerical results 

Here, we analyze CMB signals through the numerical computation of CMB power spectra, 
which are given by (4.33) and (4.34). Let us consider the case where the strength of the 
magnetic part is maximized without spoiling inflation, namely, n = 2.1. As mentioned 
in section 2, the CMB signals for this case may be enhanced to the level we can observe. 
Considering the standard single field slow-roll inflation, the slow-roll parameters are small, 
i.e., €,S <^ 1, and the tensor-to-scalar ratio denoted by r can be related with the slow-roll 
parameter as r = 16e. Neglecting the slow-roll corrections, we fix the parameters as a = 2.6 
and = /i = 3/2 and hence we have 

WjA^'^\k,,k'„k'„rj,) = 22A.575^HParkl-'k',-'-%-'-' , (4.36) 
WiA^^\kuk'^,k'^,rjj) = ^-WiA^'^\kuk'^,k'„r]i) . (4.37) 
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For explicitly obtaining the CMB power spectra, we also need concrete forms of /C^^^'s 
and K^^^^s. However, calculation of these for a = 2.6 is so complicated that we here 
substitute the expressions for a = 5/2. In appendix C, these analytical forms are provided 
and we found that IC^^^ and IC^^^ are neghgible compared with K,^^^ and K,^^^ due 
to the suppression by {—krjj)'^. In addition, we notice that in JC^^^^s and JC^^^^s, only 
}C^^^ ,)C^^^ and ifg^^ involve the terms which grow logarithmically outside the horizon 
as ln{—ujrij) with oj = ki -\- k2 + k'^. Considering the contributions on the interesting scales 
to calculate the CMB signals, namely, 10~^ ^ k x Mpc < 0.1, ln(— w?//) is just identical to 
the e-folding number and it is about —60 ~ —50. On the other hand, the other terms in 
j^TlBB jqUBB g^j^^ K2^^ have the power-law dependence on k and their coefficients are order 
of unity. Therefore, we believe that the logarithmic terms dominate over and 
j<^hBB ^ Following this concept, we use such approximate forms as 

irUBB irUBB ^ I60k2k!^ 

I\,n = — /Vq ~ — 



T^hBB ^ 4:0k!^k'^ 



(4.38) 



and neglect other /C's and K^s. 

Figure 1 shows the CMB XX (left top panel), X8 (right top one), ££ (left bottom 
one) and BB (right bottom one) power spectra generated from {TZYY) (red solid curves) 
and {hYY) (green dashed ones) respectively given by equations (4.33) and (4.34) (hereafter 
referred to as the SSa and TTa cases, respectively). For comparison, we also plot the CMB 
power spectra induced by the scalar (blue dotted curves) and tensor (magenta dot-dashed 
ones) auto-correlations, which are formulated as 

= 11 kldk^Pn{ki)T^x%^h)r^l\Ski) ' (4-39) 
C^SIa = lj kldk,Pn{k,)T^X{k,)T^^,^{k,) , (4.40) 

(hereafter referred to as the SS and TT cases, respectively). Here, Pnik) = and 
Ph{k) = 2|/ifc*p are power spectra of primordial curvature perturbations and primordial 
gravitational waves on the superhorizon scales, respectively. To obtain these spectra, we used 
the modified version of the Boltzmann Code for Anisotropics in the Microwave Background 
(CAMB) [20, 62] and the Common Mathematical Library SLATEC [63]. 

In this figure, comparing the red solid curves with the blue dotted ones, we can see 
that the overall behaviors of Cis for the SSa case are in good agreement with those for 
the SS case. This may be due to the same transfer functions of both cases. As seen in 
equations (4.33) and (4.39), both cases have completely different dependence on multipoles. 
Nevertheless, the transfer functions strongly impact on the shapes of the CMB power spectra. 
Although the tensor modes are not as clear as the scalar modes, similar behaviors can be 
observed in the TTa and TT cases. 

In contrast, there are differences of the amplitudes between the SSa and 55" cases. From 
this figure, we can see C^XiX2 ^'xiX2 1 ~ O-Ol- This value can be understood by the naive 
estimation of the magnitudes of primordial perturbations. By using equations (2.14), (2.18), 
(4.36), and observational values of the cosmological parameters [50], we gain (^^^^C^*^^) ~ 2.4x 

10-9 and (e^^^^A^ 

~ 3.3 X 10 As seen in the previous sections, since C^XiX2 1 indepen- 
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Figure 1. Absolute values of CMB power spectra of the II (left top panel), I£ (right top 
one), ££ (left bottom one) and BB (right bottom one) modes. Red solid (SSa), green dashed 
{TTa), blue dotted {SS) and magenta dot-dashed (TT) curves correspond to the spectra sourced 
from the electromagnetic-scalar bispectra (4.33), electromagnetic-tensor correlation (4.34), scalar 
auto-corrclation (4.39) and tensor auto-correlation (4.40), respectively. Here, we have adopted 
Hj = 5 X 10^^GcV,a/ = 10~^^ and r = 16e = 0.26. The other cosmological parameters have 
been fixed to the mean values reported in ref. [50]. 



dent of Wi, it can be believed that C^xfx2//CxjX2,e ~ I (C^°^C^°^> is a good approx- 

imation. Hence, we can get the consistent result. In the same manner, from equation (4.37), 
the magnitudes of primordial tensor perturbations are evaluated as {|^(^2)^(±2)^^ ^ 1.6 x 10"^'^ 

and i^S^'^^ ^^'^^^ ~ 1-3 x 10~^^. Therefore, we see that at large scales, in particular, a naive 
estimation, i.e., C^x^xIJ ^x^x^/ ~ (-^^^^^d^^^) / (-^^^^^C^^^^) ~ 0.1, is consistent with fig- 
ure 1. A fact that C^x^-^2 1/ ^'x^X2 1 about ten times larger than c'^Xixl il ^^XiX2 1 
to the difference of the metric perturbations induced by electromagnetic fields between the 
scalar and tensor modes as shown in equation (4.3). Unfortunately, the SS case dominate 
over the CMB power spectra and the SSa signals are subdominant components. Regardless 
of it, the above fact may help us to differentiate the SSa signals from the 55 ones, respec- 
tively. On the other hand, at small scales, the TTa spectra is comparable to the TT spectra. 
This enhancement can be observed in the BB mode since there are no noise of the SS spectra 
unlike in the XX^XE and EE modes. 

Furthermore, interestingly, the XX and EE modes for the SSa cases have negative signals. 
More precisely, the SSa signals are opposite in sign to the SS ones in the XX, XE and EE 
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modes. This is a consequence of the fact that the SSa signals are generated from the cross- 
bispectrum of two electromagnetic fields and one curvature perturbation and hence can take 
negative values. This may also become a clue as to the existence of the interaction between 
the scalar and vector fields (2.2). 

There may remain a concern about the contributions of CMB power spectra sourced 
from (YYYY). However, in the similar manner as the above discussion, from equation (4.3), 
these amplitudes are evaluated as 

and 

- iiT\T') I - 10- , (4.42) 

and thus we believe that total CMB signals do not drastically change. 

5 Summary and discussion 

In this paper, we studied the impacts of primordial cross-bispectra between the vector fields 
and metric perturbations on the CMB anisotropics. In the previous study [47], only the 
cross-bispectrum between one curvature perturbation and two magnetic fields is analyzed. 
However, the kinetic term of the vector field (2.2) also induces the cross-bispectra composed 
of electric fields and tensor perturbation. For the sake of completeness, we first formulated 
the primordial cross-bispectra: {TZEE) , {TZBB) , (hEE) and (hBB). Then, we found that 
{IZEE) and (hEE) {{TZBB) and (hBB)) dominate {TZBB) and {hBB) {{TZEE) and {hEE)) 
if re = —2 (re = 2). The contributions of electric and magnetic parts seem to turn over 
depending on the sign of re. We also confirmed that like the magnetic-scalar case [47], the 
cross-bispectra include the term, which expresses the logarithmic growth as ln[— (A;i + k2 + 
^s)^/] and produces significant signals. 

The anisotropic stress, which consists of the square of the electromagnetic parts of the 
vector field, acts as a source of the CMB anisotropy. Hence, {TZEE) , {TZBB) , {hEE) and 
{hBB) induce not CMB bispectra but CMB power spectra. From formulation of such CMB 
power spectra, it was confirmed that they do not have the mode-coupling components, e.g., 
the scalar-vector correlation, and are independent of the coupling constant in the kinetic 
term of the vector field (2.2). 

In numerical analysis for the case where the magnetic part is maximized, we ob- 
served that Ce's generated from the electromagnetic-scalar and electromagnetic-tensor bis- 
pectra have similar shapes to those induced by the auto-correlations of the primary non- 
electromagnetic scalar and tensor perturbations, respectively. With respect to the ampli- 
tudes, the electromagnetic-scalar spectra are about 1% of the scalar auto-correlated spectra. 
Interestingly, the former signals have opposite signs of the latter ones. In the tensor modes, 
although the signs are not flipped, the ratio between the electromagnetic-tensor and tensor 
auto-correlated spectra is improved to more than 10%. Especially, at small scales, the signals 
are enhanced and therefore the shape of the BB mode changes compared with the primary 
non-electromagnetic case. The above characteristic behaviors may inform us of the existence 
of the coupling between the scalar and vector fields (2.2). 

In this paper, we focused on the case where CMB signals are maximized, i.e., re = 2.1. 
However, their scale dependence and magnitudes are sensitive to re. This value directly 
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reflects the model of megnetogenesis; hence detailed probes with actual observational data 
remain as a future issue. 
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A Projection vectors and tensors 

Here, let us summarize the conventions of the projection tensors based on refs. [37, 57, 60, 
61]. To calculate actual CMB power spectra, we shall fix the definition of a unit vector, a 
normalized divergenceless vector and a transverse-traceless tensor. These are respectively 
given by 



'a 



m 




(A.l) 



m 




with 




(A.2) 



Here, the contraction of a™'s is easily calculated as 




(A.3) 



A divergenceless vector and a transverse-traceless tensor obey 



6W(k)ei^')(k) = 5A,-A' (for A, A' = ±1) , 
eif)(k) = fc4f)(k) = 0. 



(A.4) 



eli) (k)elf (k) = 26x,^x' (for A, A' = ±2) . 
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By a diver genceless vector and an antisymmetric tensor, a unit vector is also expressed as 

k = iVabce^a''\^)i^'\^) , (A.5) 

which is used in calculation of equation (4.12). 

Obeying these notations, we define the projection vectors and tensors of the scalar, 
vector and tensor modes, respectively, as 



of (k)EEifc, , 

o(±^)(k) = -z4^i)(k). 



(A.6) 



and 



-2/-- Y. ^i'«<-"rvM„„„, 

Mnianif, 



2 1 1 



= ±2V34\-^ E Ti>^2M(k)a^ar ( 

Mniami, 

'211 



2 1 1\ (A.7) 

M nia mi, 



^ \ M rua ruf, 



where I2I1 ""^ = y ^ is given by equation (4.11). These decompose arbitrary physical vector 
and tensor such as the velocity, the metric and the anisotropic stress into the scalar, vector 
and tensor components: 

coaik) = u;(0)(k)OW(k) + E ^(')(k)Oi^)(k) , 

A=±l 

Xabik) = -^Xisoimb + X^°Hk)O^Jik) (A.8) 

A=±l A=±2 

From equation (A. 4), we have the inverse formulae as 

a;W(k) = -Of(k)a;,(k) , 
u;(±i)(k) = -0(Ti)(k)o.,(k) , 

x(°Hk) = ^OS)(k)Xa.(k), ^^g^ 
X(±^)(k) = i0if)(k)xa.(k), 

x(^^Hk) = ioif)(k)x..(k). 
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B Time-dependent mode function 



Here, we explain properties of the time-dependent part of the mode functions, Um{x), given 
by equation (2.12). 

For specific m, we have the analytic formulae as 



^^5/2(2;) 



^(-3 + 3ix + x2)e^^ 



U3/2{x) = (1 -«x)e*^ , 



-3/2 c 



The derivative with respect to x is given by the recurrence formula as 

dum{x) 2m 



dx 



(x) - Um+l{x)] . 



In particular, we have 



^U5/2(x) = -^x(-l + ix)e* 



d 



-n_3/2(x) = 3(^ + ^-^ 



3i 



x^ 



At the limit: x — )■ 0, these asymptotically behave as 

Um(x 0) 
'dum{x) 



dx 



1 (m > 0) 

(m = 0) , 

K r(m)r(m+i) ^ ym <.u) 

2i~^'" ^(»-cotm^) ^2m-l < 1) 



r(m)2 



C Specific expressions of K^^^ and /f^^^ 

C.l u = fi = 3/2 and a = 5/2 

The time-integlas in equations (3.7) and (3.14) are 



j ^ driT^ '^U3/2iTi) (^^5/2 (^2n)^ 



= i4)^)+0(-A:i7?) , 

X) 

(iTirf^U3/2(Ti)u5/2(x2Ti)n5/2(x3ri) 



(B.l) 
(B.2) 
(B.3) 

(B.4) 

(B.5) 
(B.6) 



(B.7) 
(B.8) 



(C.l) 



(C.2) 
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3(-fcir?)3+ 6(-A;i7?)A;2 +6 +'"^5/2+^^ ^^^^ 



where 



j: 



(la) 
5/2 



J, 



(Ife) 



J. 



J. 



[2b) 
5/2 



j: 



(3a) 



J, 



driT^ ^ ('^^^^^/^(n)^ '"5/2(a;2n)^i5/2(2;3Tl) 

drirf^U3/2(Ti)u5/2(x2ri) (^Ti-^n5/2(x3ri) 



(iTirf^U3/2(Ti) (^'^1^^5/2(2;2Tl)^ ^i5/2(2;3n) 



(C.3) 



(C.4) 



(C.5) 



^2 



3(-fcir?)A:2 ^ 5/2 



^21.2 



5/2 9fcfa;3 
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5/2 



^2^2 

2 3 +^(^k2 + k:i) + 2k2k^] 



(C.6) 

(C.7) 
(C.8) 



-3 [7 + ln(-w?7)] 



+ 



[ArpiCA;! + a;) - 3tjA:2fc3(a;^ - kiuj + fci) + ?,kiu?{2ki - uo) + uj^] ,(C.9) 



[7 + 



+ 



1 



[9a;^ - 3w^(A;| + Ak2hi + + 3ujkik2k3 - 2A:p|] , 



5/2 9A;3a;3 
(3b) 



[3uj^ - 3{k2 + /C3)w3 + (A:2 + 3fe 



2\ 2 



+/c2(/ci + k2k3 + 3/c|)a; + 2k'^k3{k2 + A;3)] , 
'5/2 = 4%^ik2^k3) . 
The other parts are expanded as 



(C.IO) 



(C.ll) 
(C.12) 



^3/2 (- ^1^)^5/2 (-^"2?/)n5/2 (-A;3??) 



= 1 + 1 (3A:2 + kj + kl) rj' - '-{-k^r^f + 0{{-kriif) 

^/2 (-^1^)^5/2 (-^2?/X/2 (-^3??) 

= ^ [1 + + kl + kiw + 1(14 + kl + ki)v' 



(C.13) 



+ Oii-k,vf) ■ (C.14) 



Picking up the leading-order terms in the combinations of these equations, we can reach the 
final expressions of i^^^^'s and K'^^^^s at the end of inflation on the superhorizon scales 
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{-krji < 1): 
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C.2 zy = ^ = 3/2 and a 

grals and the 

(iTiTiM3/2(ri) 



-3/2 



The time integrals and the product of the mode functions are resectively expanded 
4 , ^ f d , \ ( d 



as 



^^^--3/2(x2ri) 



nk\ 



27kf{3kl + k1 + k1) 27TTkf 
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81 
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In the same manner as the above case, we can gain 
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